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Abstract 

We consider the Nambu and Hamiltonian representations of Rayleigh- 
Benard convection with a nonhnear thermal heating effect proportional to the 
i-^'. Eckert number (Ec). The model we use is an extension of the classical Lorenz- 

63 model with 4 kinematic and 6 thermal degrees of freedom. The conservative 
parts of the dynamical equations which include all nonlinearities satisfy Liou- 
ville's theorem and permit a conserved Hamiltonian H for arbitrary Ec. For 
Ec = two independent conserved Casimir functions exist, one of these is as- 
sociated with unavailable potential energy and is also present in the Lorenz-63 
truncation. This Casimir C is used to construct a Nambu representation of 
(— I ■ the conserved part of the dynamical system. The thermal heating effect can be 

fS f represented either by a second canonical Hamiltonian or as a gradient (metric) 

system using the time derivative C of the Casimir. The results demonstrate the 
impact of viscous heating in the total energy budget and in the Lorenz energy 
cycle for kinetic and available potential energy. 

Keywords: Rayleigh-Benard convection, Nambu mechanics, Lorenz equations. 
Viscous heating 



1. Introduction 

The Nambu representation is an extension of Hamiltonian dynamics with 
the aim to formulate dynamical equations satisfying the Liouville theorem (20j . 
In state spaces with dimension above two further conserved quantities besides 
the Hamiltonian can be used. Thus Casimir functions can be considered as 
additional Hamiltonians. This concept is useful in dynamical systems when 
a physical interpretation of total energy is not available. Nambu systems in- 
corporate noncanonical Hamiltonian systems with the conservation laws being 
Casimirs of the Poisson tensor. Algebraic properties are formulated in terms of 
anti-symmetric and cyclic extensions of Poisson brackets. The development of 
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Hamiltonian low-order models with additional conservation laws has been also 
put forward in terms of gyrostat concept 0, 0]. The Nambu representations 
of hydrodynamics (2]| has shown up to be useful in the design of numerical 
algorithms in geophysical fluid dynamics since they deliver a general tool to 
preserve conservation laws 2^ 2^ ■ The Casimirs of ideal hydrodynamics have 
been explicitly derived in the noncanonical Hamiltonian representation ■ The 
application of Nambu mechanics in meteorology was put forward in recent years 



The continuous Rayleigh-Benard convection can be recast in a Nambu sys- 
tem using total energy and a Casimir that is considered resulting from Kelvin's 
circulation theorem [2| or as 2D analogue of helicity 27 1. Viscous dissipation 
is present in a symmetric bracket using the difference between total energy and 
unavailable potential energy. Symmetric brackets have been proposed as a nat- 
ural extension of symplectic brackets and the combined system has been coined 
'metriplectic' [l^ (combining metric and symplectic systems). 

Saltzman derived the Fourier mode truncations of the 2D Rayleigh-Benard 
convection in the Boussinesq approximation [30| . The three- mode truncation 



of these equations in the classical Lorenz-63 model [12| includes a mode for 
the stream-function, a mode for temperature perturbation, and a mode for the 
vertical stratification. Lucarini and Fraedrich [3] included heating by viscous 
dissipation characterized by the Eckert number Ec in a low-order model with- 
out shear flow. In order to incorporate this effect, the Lorenz-63 model has 
to be extended to 10 degrees of freedom including four stream function, four 
temperature perturbation, and two vertical stratification modes. 

Some of the results presented in suggest that the consideration of the 
Ec terms is relevant. First, it was found that when Ec=0 the system obeys a 
symmetry, which, by Noether's Theorem, corresponds to a conserved quantity. 
Such a symmetry is broken when positive values of Ec are considered, with the 
result that the phases of the temperature and streamfunctions are efficiently 
mixed and degenaracies are removed in the system, so that only one vanishing 
Lyapunov exponent is found. In flT] , it was proposed that neglecting the heating 
effect due to viscous dissipation could be responsible for the observed imbalance 
of the global energy balance of the state-of-the-art global climate models. 

The dynamical system given by the Lorenz-63 equations can be decomposed 
in a conservative part which conserves volume in phase space and a divergent 
part related to forcing and dissipation. The conservative part has two con- 
served quantities, a Hamiltonian H and a Casimir C, which can be used to 
determine the dynamics in a three-dimensional Nambu system ^ m so that 



{X, Y, Z) = VC X ViJ + V$, where we have used the classical notation. Forc- 
ing and dissipation is included in the gradient of $. Thus Nambu dynamics 
provides a clear geometric view of the nonlinear dynamics and yields hints to 
the geometry of the strange attractor in the forced and dissipative equations 



22, 1^. The Nambu structure of a six mode truncation of Saltzman's [Sil con- 



vection model has been presented recently [3,] . The derived model has 5 modes 
in common with the model studied here. 

The aims of this publication are (i) to present the Nambu representation of 
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the conservative part of the 10 degrees model of Rayleigh-Benard convection 
presented in [l^ and (ii) to determine the algebraic properties of heating by 
viscous dissipation. In Section [5] the truncated dynamical system is presented. 
The conservation laws are derived in Section [3] with a physical interpretation in 
Section The dynamics is represented in terms of a Nambu system using two 
conservation laws (Section [5]) from which a noncanonical Hamiltonian structure 
can be derived (Section [6]). The viscous heating effect is shown to possess a 
symplectic as well as a metric structure (Section [7]) . Results are summarized 
and discussed in Section [51 

2. Truncated equations 

The physical setting is the two-dimensional Rayleigh-Benard convection 
driven by a vertical temperature contrast AT across the height H. Dynamics in 
the vertical plane (a;, z) follows the Boussinesq approximation for incompress- 
ible velocity u = —d^ip, v = d^-ip^ with the stream-function ip, and the linearized 
equation of state p = poi^ — oT) with the thermal expansion coefficient a. 
Temperature is decomposed according to T = Tq — z/ST / H + in a mean tem- 
perature To, a linear temperature gradient and the perturbation 6. Boundary 
conditions at z = and z = H along the x-axis demand vanishing vertical 
velocity, -0 = 0, vorticity, V^?A — 0, and temperature perturbation, 6 — Q. 

The main extension of the Boussinesq equations suggested by [l^ is the 
inclusion of viscous heating in the temperature balance 

5f V^V- -I- J(V', V^V') = gad^e + vV^iP (1) 

dtO + j(V', e) - ^a.v- + kv^o + vs.jV (2) 

here v is viscosity, k heat conductivity, is heat capacity, and g is gravity 
acceleration. Non-dimensional equations are obtained hy x — Hx' , t = H^t' /k, 
tP = kip' and e = kuO'/gaH^. 

The impact of heating by viscous dissipation in the Fourier decomposition 
requires at least two modes for the stream-function and the tempera- 

ture perturbation 9n,m, furthermore, two modes for the perturbation of the 
stratification need to be included, 0o.2 and 6*0,4 jlH. The Fourier modes are 
associated with the wave vectors (Tram, tto). Shear flow is excluded here by 
^'o.i = 0- Real and imaginary parts of the stream-function are split according 
to ^1.1 = a{Xi -\- 1X2), 5*2.2 = Ci{Ai + 4^2), and the temperature perturbation 
is 6ii,i = p{Yi + iY2), = P{Bi + 1B2), Baa = i^Zi, and 604 = ^7^2- 

Heating by viscosity is proportional to the Eckert number, Ec = 
k^ / (CyATH^) and leads to new nonlinear terms in the temperature equations 
of the Lorenz system. The present analysis concentrates on the conservative 
parts of the dynamical equations with zero divergence in phase space. The con- 
servative parts are obtained by an operator splitting through the neglect of the 
diagonal linear forcing terms (see the general decomposition (4) in |12i). 
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The conservative parts of the truncated equations retain ah nonhnearities 



Xl - aY2 (3) 

X2 = -aVi (4) 

Yi = X2Zi~rX2 + 5/{V2TT)aEc{XiAi+X2A2) (5) 

Y2 = -XiZi+rXi+5/{V27r)aEc{XiA2- X2A1) (6) 

Ai = a/2B2 (7) 

A2 = -a/2Bi (8) 

Bl = AA2Z2-2rA2-l/{2V2Ti)aEc(Xl- xl) (9) 

B2 = -AAiZ2 + 2rAi-l/{V2TT)aEcXiX2 (10) 

Zl = Xiy2-^2il (11) 

Z2 = 4:AiB2 - 4:A2Bi (12) 



The non-conservative terms not considered here can be represented as a hnear 
gradient, — Va;$, of the potential 

$ = a{Xf + Xi + AAl + AAl) + Y^ + + + B| + feZj^ + 46Z| (13) 

with respect to the 10 variables ([3]) - p^ . see [l^l- The notation follows the 
standard Lorenz-63 model with the Prandtl number <j = v /k and the relative 
Rayleigh number r = R/Rc, with i? = gaH^AT/kiy, Rc = 277r^/4, and a 
geometric factor b; for further details see [l^, [l3|- The symmetries and the 
instability of these equations have been intensively studied in [l6| . 



3. Conservation laws 

Due to the absence of linear forcing and dissipative terms in ([3]) - (IT2|l the 
phase space divergence vanishes and the equations satisfy the Liouville theorem. 
The equations have several time independent conservation laws (determined by 
the Reduce package CONLAWl 

The first conservation law which is preserved for arbitrary Eckert number Ec 
represents the sum of kinetic and potential energy and is denoted as Hamiltonian 

H = ^ {Xf + Xl + 4Al + 4Al - 2aZi - aZa) (14) 

Further conservation laws are found for Ec = 0; The first is quadratic in the 
temperature perturbations 

C = i (Yj^ + Yi + Bl + B'I + zl + zl - 2rZi - rZz) (15) 

and the second is based on a coupling between stream-function and temperature 
perturbation 

C =]^{XiYi+X2Y2 + AiBi+A2B2) (16) 
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The conserved quantities C and C are denoted as Casimir functions according 
to there specific interpretation in Hamiltonian dynamics (see Section [6]). 

The Haniiltonians and the Casimirs for the two wave numbers considered are 
conserved individually (pairs of variables X,Y for {n,ni) = (1, 1), and (-A i?) 
for (2, 2)). The sum of H and C has been analysed in a related truncation [6ll34| 
(denoted as "Q"). 

The time derivative of the Casimir C is linear in Ec 

+ 10Y2{A2Xi - A1X2) - Bi{X^ ^ X^) - 2B2X1X2] (17) 

The time derivative of the second Casimir C is determined by the iJc-dependent 
dynamics of Bi and B2 only 

^ = ^ (^1 ix! -Xl) + 2A2X,X2) (18) 

The time derivative of C in ([TT]) will be used to formulate the i?c-dependent 
parts in terms of a gradient system (see the forthcoming Section [7^ . 

4. Physical interpretation of the conservation laws 

The total energy is the sum kinetic and potential energy 

H = KE + PE (19) 
which are given as integrals in the vertical plane 

KE^^J {VijfdV (20) 
PE = -ag I z9dV (21) 



(see for example [5|). 

In the present truncation the kinetic energy KE is given by 

KE=^ {Xf + Xl + + AAl) (22) 

and total potential energy is 

PE = -^{2Zi + Z2) (23) 

The available potential energy characterizes the amount of potential energy 



which can be converted to kinetic energy [111 [3 



APE^^^Ie'dV (24) 
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The truncated version is 



APE = - ^ {Y^ + Yi + Bf + Bl + Zl + Zl) (25) 

Note that APE < since any perturbation reduces the energy transformable 
to kinetic energy. Total energy and the Casimir C are related to KE and APE 

by 

H = KE + APE + -C (26) 

r 

Therefore, the Casimir contribution can be identified as unavailable potential 
energy Q 

-C ^PE~ APE (27) 
r 

The balance between kinetic and available potential energy is given by the 
Lorenz energy cycle 

j^KE = G {X1Y2 - X2Y1 + 2A1B2 - 2A2B1) (28) 

^ape^-IkE^^^ (29) 
dt dt r dt ^ ' 

For finite Eckert number {Ec > 0) viscous heating modifies available potential 
energy (see dC/dt in (|f 7p ). The balance between KE and APE is closed, 
dKE/dt + dAPE/dt = 0, in the absence of viscous heating. The sum E of 
kinetic and available potential energy is the available energy 

E = KE + APE (30) 

and will be used represent the viscous heating term (see ([38| below). 
The Casimir functions (IT51 ITC)) correspond to the 

Continuous versions of the Casimir functions are 0, [l^ Ci = J g(T — z)dV 
and C2 = J Qh{T — z)dy , where T is the scaled deviation of the temperature 
from a linear profile z, C, is vorticity; the integrals are in the vertical plane. Ci 
is based on the preservation of T — z contours and corresponds to C in (fT5|) 
for the function g{x) = jl. C2 depends on Kelvin's circulation theorem and 
corresponds to C in (ITB)) . In the continuos Nambu representation of ^ the 
Casimir C2 is used together with the Hamiltonian H ^ J (1/2(V'0)^ — Tz)dV, 
while [i^] compare different Nambu representations using H, Ci (for h{x) = x), 
enstrophy and total buoyancy (Ci with g{x) = x). In the present analysis the 
first Casimir C in pS]) . i.e. the unavailable potential energy is used in order to 
continue previous work on the Nambu representation of the Lorenz-63 model 
(2^ . A further motivation for using this Casimir is its role in the viscous heating 
(see Section [7T2|) . 
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5. Nambu structure 



Nevir and Blender 22| developed a Nambu representation of the non- 
dissipative parts of the classical Lorenz-63 [l2| model for Xi,Y2, and Zi (see 
also jlj, ^] and ^25;] for a classification). With the conservation laws for H 
and C the Lorenz equations are in terms of a; = (-'^i, Zi) 

= -EijkCjHk, i,j,fc = 1,2,3 (31) 



with the Levi-Civita symbol e^fc (note that the definitions in 22| are inter- 
changed). The second Casimir C does not occur in the classical Lorenz-63 
truncation since Yi and X2 are not considered. 

Rewriting the dynamical system ([3]) - (fT2)) for the 10-component variable 
vector 

x^{X,,Y2,Z,,X2,Y,,A,,B2,Z2,A2,B,) (32) 
the Nambu equations for the extended Lorenz system is 

dx 

— -AA(V,C,V,iJ) (33) 

with the Nambu tensor JV and the gradient Vx- 

The non- vanishing elements of the Nambu tensor are 

^1,2,3 ~ — £l23: -^3,4,5 = £345; -^6,7,8 = —£678, -^8,9,10 = £8910 (34) 

with the standard definition of the Levi Civita symbol e. Thus the Nambu 
tensor is cyclic and anti-symmetric. Furthermore, it satisfies the generalized 
Jacobi identity |33l, |31| 

N.jkNipq + N,jqNikp + N.jpNiqk = (35) 

in the form required for constant Nambu tensors. 

Based on ([55]) a Nambu bracket for arbitrary phase space functions F{x) is 
written 

— = {F,C,H} (36) 

which is cyclic {F,C,H} — {C,H,F} = {H,F,C}. This cyclicity is the main 
ingredient in applications of the Nambu representation in numerical models [i^ . 
This bracket vanishes for the second Casimir (|16|) . {C,C,H} = 0. 

The Nambu representation of the conservative terms can be presented in 



7 



terms of brackets 
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where the bracket {,} is the Jacobian of C and H with respect to a and 6, 
{a, 6} = dC/da — dH/db. The three brackets in the upper left represent the 
Nambu dynamics of the Lorenz-63 ([5T|) model [221 . This structure demonstrates 
the decoupling of the two wave number modes (to, n) = (1, 1) and (2, 2) in the 
five upper and the five lower rows. Furthermore, the coupled Fourier modes are 
separated in the columns. Coupling between the wave number modes is given 
by the feedback represented in the Eckert terms which is not included here (see 
Section [T]). 

Instead of the Hamiltonian the dynamical equations can be rewritten using 
the available energy E = KE + APE in the Nambu bracket (p6| 



since E deviates from H only by the Casimir 



(38) 



H = E+-C 
r 



(39) 



which is annihilated in the bracket. 



6. Hamiltonian structure 

The Hamiltonian structure of the dynamical system ([3]) - ((T^ is a straight- 
forward evaluation of the Nambu system (1551) (also for Ec = 0) 

j^x=VV^H =N{V.,C,V^H) (40) 

with the anti-symmetric Poisson tensor 

7'(.)=AA(V,C,.) (41) 
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The elements of the Poisson tensor are exphcitly given by derivatives of the 
Casimir with respect to the variables in (|32]) . for example Ci — dC/dXi, etc. 
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(42) 

The Poisson tensor V defines a Poisson bracket for the dynamics of arbitrary 
phase space functions F{x) 

'i-{F,H}. (43) 
The Casimir C is conserved as the nuUspace of the Poisson tensor 

V\/^C = (44) 

Thus, in Hamiltonian dynamics the conservation of Casimirs is an inherent 
property hidden in the Poisson tensor. These conservation laws lead to foliations 
in phase space and are interpreted as a kinematic property of the dynamical 
system, whereas the conservation of the Hamiltonian is considered as a dynamic 
property. 



7. Viscous heating by Eckert terms 

An important property of the i?c-dependent terms in the equations for 
Yi,Y2,Bi, and B2 ®, ®, and ^ is that they satisfy the Liouville theo- 
rem. Thus a representations in terms of a stream-function is possible (hence as 
Hamiltonian systems). To obtain this we consider the equations for Yi,l2 and 
Bi, B2 in two complex plains given hy Y — Yi + iY2 and B = Bi + iB2. The 
dynamic equations are determined by the gradient of the complex function 

f{Y,B) = -^^{lOX* AY* - XX B*) (45) 
2v27r 

with X = Xi+iX2 and A = Ai+iA2. The Cauchy-Riemann equations yield the 
dynamics for Y and B in terms of a stream-function (Im /) and the gradient of 
potential (Re /). It turns out that the potential is given by the time derivative 
of the Casimir C (unavailable potential energy) which leads to a metric system 
(loj . The stream- function is considered as a second Hamiltonian H (not to be 
confused with ([M]) ). thus the complex function is / = (7 + iH. In the following 
the Hamiltonian and the gradient representation are considered. 
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1.1. Hamiltonian system 

The Hamiltonian H — Imf ([jS)) can be constructed by the Ec-dependent 
forcing terms for the set of variables yj — {Y2,Yi,B2, Bi} in ([5^ 

J^y, - G[y,] (46) 

using that these terms G depend only on Xi, X2, Ai, A2. Therefore, H is for- 
mally 

H = Y2G[Yi] - YiG[Y2] + B2G[Bi\ - BiG[B2] (47) 
which can be rearranged in the compact representation 

~ n d Yi r, d Bi , . 

H^Yi -+Bl (48) 

^ dtY2 ^ dt B2 ^ ' 

determined by the phases of 6*1.1 — PiYi + iY2) and 6*2.2 — P{Bi + iB2). In 
terms of the dynamical variables, H has the explicit form 

H = ^^[10^2(^1X1+^2X2) 
2v 27r 

- 10^1(^2X1-^11X2)- 52(^2 -X|) + 251X1X2] (49) 

Thus the viscous heating effect can be represented by a second Hamiltonian 
part 

_d 

dV 

where Vj, is the gradient with respect to y. The Poisson tensor is canonical 
with the non-zero elements 



-y = VEcVyH (50) 



Vec{Y,,Y2) = -1, pEciY2,Y,)^l 

Vec{Bi,B2) = -1, Vec{B2,B^)^1 (51) 

The dynamics ([50)1 has a canonical Poisson bracket 

{F, H}ec - VxFVecVxH (52) 

which conserves the Hamiltonian (ITil) . {H, H}ec ~ 0. Thus H can formally be 
considered as a Casimir of the tensor Vec, VecH = 0. 

7.2. Gradient System 

The i?c-dependent terms in (jSj - ([T2|) are the gradient of a potential given 
by the real part Re/ (|45p . This potential can be identified as time derivative 
of the Casimir ()17|) . Hence the present low order system shows a relationship 
between the Hamiltonian part ([^^ and the viscous heating mediated by the 
Casimir ([15]). The potential can be constructed from the Ec-dependent forcing 
terms (|46| according to 

dG 

— = riG[yi] + y2G[r2] + b^g[b{\ + B2G[B2\ (53) 
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Embedding this in the dynamics of x the gradient representation is 

with the symmetric operator M. = diag{0, 1, 0, 0, 1, 0, 1, 0, 0, 1) which selects the 
variables Y2,Yi, Bi, B2 in the complete 10-component vector x (|32p . Thus for 
an arbitrary phase space function F the contribution of viscous heating (j54p can 
be included in a symmetric bracket (also denoted as metric bracket in contrast 
to symplectic Poisson brackets). 

(f.;|,.V,FA.V.f (55) 

Due to the Cauchy-Rienimann equations the H manifolds are perpendicular to 
C and H is preserved in this bracket, {H,C) = 0. With (1551) the dynamics of 
the full system can be completely understood in terms of the available energy 

since 

dC r dE 

-r ^ ;- (57) 

dt a dt ^ ^ 

represents the change in available energy. 

This symmetric bracket does not represent physical dissipation associated 
with the convergence of trajectories in phase space. The reason is that the dis- 
sipated energy is re- introduced in the temperature equation which is represented 
by dynamical degrees of freedom. 

A main property of this symmetric bracket is that it generates a positive 
change of the Casimir tendency 

(C, C) = (100(A? + Al){Xf +xl) + {Xi + Xif) > (58) 

note that this is proportional to Ec^. Since the Nambu bracket representing 
the conservative part does not vanish {C, C, E} 7^ 0, the time derivative of the 
Casimir tendency in the full equations has no definite sign 

^^C = {C,C,E} + {C,C} (59) 

Due to this property, C cannot be considered as the analogue to entropy (see 
the concept in [9]). 

Since the Hamiltonian is conserved for arbitrary Eckert number the sym- 
metric bracket vanishes for H 

— ^{H,C}=0 (60) 

Hence the Hamiltonian can be considered as a nuUvector of the symmetric 
bracket (see the construction method (l9j). 
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The role of the Casimir is twofold: (i) As a conserved quantity it deter- 
mines the non- viscous dynamics in the Nambu representation, and (ii) the time 
derivative determines the heating effect in a gradient system. 



8. Summary and Conclusions 

This paper discusses the impact of viscous heating on the dynamic properties 
of convection. The model is given by a truncation of the Boussinesq equations 
(soj with an additional heating caused by viscous dissipation [l^ . Wave numbers 
selected are n = m = 1 and n = m = 2 for the stream function and the 
temperature perturbation, where additional vertical wave numbers m = 2 and 
m — A are included. Vertical shear is not incorporated. Linear forcing and 
friction are removed in the dynamical equations to obtain a divergence free 
dynamical system satisfying Liouville's theorem. 

The truncated systen shows one conservation law for all Ec which is given 
by the sum of kinetic and potential energy and identified as the Hamiltonian of 
the dynamical system. There are two quadratic conservation laws for vanishing 
Eckert number, Ec = 0, a first, C, related to unavailable potential energy (UPE, 
the residuum between total energy and available potential energy) and a second, 
C, related to the coupling of stream-function and temperature perturbations. 
Both conservation laws are Casimirs of the Hamiltonian representation. The 
viscous heating changes both Casimir functions. An available potential energy 
APE is identified in the truncated system which is convertible to kinetic energy. 
The Lorenz energy cycle for kinetic and available potential energy has a source 
term linear in Ec. 

A Nambu representation is defined using the Hamiltonian H and the Casimir 
C to determine the time derivative of an arbitrary state space function by a 
Nambu bracket, F = {F, C, i?}. The Nambu bracket is cyclic and satisfies 
a generalized Jacobi identity. In a Nambu representation the Casimirs can 
be considered as additional Hamiltonians. The Nambu representation can be 
directly transformed into noncanonical Hamiltonian dynamics by evaluating the 
Casimir dependency in the Nambu bracket, {F, C, H} = {F, i?}p. 

Casimirs and the Hamiltonian are excellent observables for the study of 
the response of dynamical system to small perturbations far from equilibrium 
13 ! 14 1 . Even when these quantities are not conserved in the full nonlinear 



equations, they provide an intrinsic skeleton characterizing the geometry of the 
phase space trajectories. Promising results have been obtained recently for 
perturbations of the Lorenz 96 model considering total energy and momentum 
analogues [l8l |. 

The final result of the present publication is that the dynamics Q - (|12p 
can be decomposed in a Nambu representation and a metric representation of 
the viscous heating which is proportional to Ec 

F = {F,C,H} + {F,C) (61) 

The Casimir C appears in the Nambu as well as in the metric term (via the 
time derivative). Furthermore, the viscous heating bracket is equivalent to a 
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canonical Hamiltonian bracket, {F, C) — {F, H}ec, due to the Cauchy-Riemann 
equations. 

The present low order model revels specific properties of the viscous heating 
which is usually neglected in models: 

1. Viscous heating does not perturb conservation of total energy (the Hamil- 
tonian). This result may have implications for considering this effect in 
complex models. 

2. The viscous heating terms are divergence free in phase space, hence the 
full dynamical system satisfies the Liouville theorem. 

3. The Lorenz energy cycle for kinetic and available energy has a source 
term given by viscous heating. This effect alters convective processes by 
modifying available potential energy. 

4. Viscous heating determines the time derivative of the Casimir in the 
Nambu representation. Nevertheless, its role in the Nambu and the Hamil- 
tonian representation is unaltered. 

5. The viscous heating terms can be represented by a symmetric bracket with 
a positive change of the Casimir time derivative, (C, C) > 0. However, 
the interpretation of C as an entropy according to [9| is not possible since 
the Poisson bracket does not vanish, {C, H}p ^ 0. 

This study shows that a correct energy recycling budget leads to conserva- 
tion of total energy and retains the Hamiltonian and Nambu structures. The 
single impact is a change of the Casimir function, without altering it's role in the 
Nambu representation (such conservation laws have been termed constitutive in 
Nambu dynamics [1^ 113]). The role of the Casimir in the recycling process 
might be a particularity of the present model, but the feedback in the Lorenz 
energy cycle and the recreation of kinetic energy could shed new light on the 
relationships with entropy [^, [l9| . The Ec-term mixes available potential energy 
(APE) and unavailable potential energy. This implies that an im p ro p er treat- 
ment of the Ec-term can have a large impact on the energy cycle Il l |.|l7 | and bias 
considerably the thermodynamic efficiency of a geophysical fluid [ITITtI H, [l^ . 
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